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Abstract. Graph dynamical systems (GDSs) can be used to describe a wide range of dis- 
tributed, nonlinear phenomena. In this paper we characterize cycle equivalence of a class 
of finite GDSs called sequential dynamical systems (SDSs). In general, two finite GDSs 
are cycle equivalent if their periodic orbits are isomorphic as directed graphs. Sequential 
dynamical systems may be thought of as generalized cellular automata, and use an update 
order to construct the dynamical system map. The main result of this paper is a charac- 
terization of cycle equivalence in terms of shifts and reflections of the SDS update order. 
We construct two graphs C(Y) and D(Y) whose components describe update orders that 
give rise to cycle equivalent SDSs. The number of components in C(Y) and D(Y) is an 
upper bound for the number of cycle equivalence classes one can obtain, and we enumerate 
these quantities through a recursion relation for several graph classes. The components 
of these graphs encode dynamical neutrality, the component sizes represent periodic orbit 
structural stability, and the number of components can be viewed as a system complexity 
measure. 



1. Introduction 

Sequential dynamical systems (SDSs) were introduced in 0, [Tlj]. These are dynamical 
systems constructed from (i) a finite undirected graph Y where each vertex has a state, 
(ii) a sequence of vertex functions, and (Hi) a word w over the vertex set of Y. The SDS 
map is constructed as the composition of the functions in the order specified by w. As 
such, they represent a useful framework for describing distributed phenomena with causal 
interactions. This paper is about cycle equivalence of finite graph dynamical systems, 
which we study in the context of SDSs. Two SDSs are cycle equivalent if their periodic 
orbits are isomorphic as directed graphs. We will study how the update order affects the 
structure of the periodic orbits, and thus the long-term behavior of the system. As an 
example, we show the surprising result that if the GDS base graph is a tree then there is 
only one possible periodic orbit configuration, and this holds for any fixed choice of vertex 
functions. SDS and SDS-like algorithms occur in many application areas such as [1, E2|, 
and our results will provide a behavioral complexity measure for these. 
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This paper is organized as follows. In Section [2] we describe SDS related terminology and 
relevant background results from ll|, ll3j| on functional and dynamical equivalence of SDSs. 
In Section[3l we prove one of the main results of this paper: any two SDSs where the update 
orders differ by a cyclic shift are cycle equivalent, and this holds for any choice of vertex 
functions. Additionally, when the vertex states are taken from F2 = {0, 1}, which is the 
standard choice in most studies of cellular automata, then reflections of the update order 
also encode cycle equivalent SDSs. We also show how shifts and reflections of update orders 
have a natural interpretation in terms of source-to-sink operations on acyclic orientations of 
the GDS graph. In Section H] we introduce the graphs C(Y) and D(Y) which form the basis 
for our analysis and characterization of cycle equivalence over general graphs. These graphs 
are examples of neutral networks, and we characterize some of their structural properties. 
We study the functions k(Y) and S(Y), which count the connected components of C(Y) 
and D(Y), respectively. We show how S(Y) is given in terms of k(Y) and give several 
results for the computation of k(Y) with implications to dynamics. These functions can be 
regarded as a measure for system complexity since they are upper bounds for the number 
of SDS maps up to cycle equivalence achievable through variations of the update order. As 
a computational example we demonstrate how k(Y) increases from 0(n) for radius-1 rules 
(the elementary cellular automaton rules) to 0(n-2 n ) for radius-2 rules. We also show how 
the presence of symmetries in the base graph may allow for significantly improved bounds 
in certain cases. In the summary section we show how cycle equivalence of SDSs is closely 
related to Coxeter theory. Some of the results that we prove in this paper have a natural 
analog when translated into the language of Coxeter groups. This opens the door to use 
the rich mathematical tools and results of Coxeter theory to study sequential dynamical 
systems, something that has never been done before. 



2. Background and Definitions 

Let Y be a finite undirected graph with vertex set v[Y] = {1, . . . , n}, and edge set e[Y]. 
Since most graphs in this paper are finite and undirected we simply refer to this class 
of graphs as "graphs" and specify if needed. The 1 -neighborhood of vertex v in Y is 
Bi(v;Y) = {V G v[y] I {v,v'} G e[Y]} U {v}, and the ordered 1 -neighborhood n[v] of v is 
the sequence of vertices from Bi(v; Y) ordered in increasing order. The degree of vertex v 
is written d(v). Each vertex v is assigned a state y v G K where K is a finite set. In the 
following y v is called a vertex state and the n-tuple y = (yi, . . . , y n ) is a system stated We 
write 

(2-1) y[v] = (y n [v](l), ■ ■ ■ , Un[v](d(v)+l)) , 

for the restriction of the system state to the vertices in n[v], and let y'[v] denote the 
same tuple but with the vertex state y v omitted. The finite field with q = p k elements is 
denoted ¥ g . 



In the context of, e.g. cellular automata a system state is frequently called a configuration. 
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Let fy := (/j)iev[y] be a sequence of vertex functions ff. K d ^ +l — > K, and define the 
sequence of Y -local functions #y := (Fj) iev [y] with F^: K n — > K n by 

(2.2) Fi(yx, ■■■ ,y n ) = (y 1: yi _ 1: fi(y\i\), y i+i , . . . , y n ) . 

It is clear that fy completely determines 3V, and vice- versa. However, there are settings 
when it is easier to speak of one rather than the other. 

Let Wy denote the set of words over v[Y]H Words are written as w = (wi, w 2 , • • • ,w m ), 
w = W1W2 ■ ■ ■ w m , w = (w(l), w(2), . . . , w(m)), etc. The subset of Wy where each element 
of v[y] occurs exactly once is denoted Sy. The elements of Sy may thus be thought of as 
permutations of v[Y]. The symmetric group S n acts on system states by 

(2.3) 7 • {yi, ...,y n ) = (y 7 -i(i), . . . ,y 7 -i (n) ) . 

Definition 2.1 (Sequential dynamical system). A sequential dynamical system (SDS) is a 
triple (Y, #y, w) where Y is a graph, $ Y = (Fj)j £v [y] is a sequence of Y-local functions, and 
w = (wi, . . . ,w m ) G Wy. The associated SDS map [3y,w]: K n — ► K n is the function 
composition 

(2.4) [3y, w] = F tom o F^., o-oF ffl2 oF M1 . 

The graph K of an SDS is called the base graph, and when w G S'y , the SDS is referred to as a 
permutation SDS. A sequence of Y-local functions is Aut(Y) -invariant if 7oF„ = F^^cy 
for all u G v[y] and all 7 G Aut(y). Here, the composition of a function — * K n with a 
permutation of K is interpreted as in (I2.3p . The corresponding sequence of vertex functions 
fy is Aut(Y)-invariant if 3v is Aut(y)-invariant. The phase space of the map cf>: K n — ► K n 
is the directed graph T(0) with vertex set K n and edge set {(y,0(y)) | y G K n }. The 
following example illustrates these concepts. 

Example 2.2 (Asynchronous Elementary Cellular Automaton rule # 1.). Let Y = G1X4 
which is the graph with vertex set {1, 2, 3, 4} and edges all {i, i + 1} with indices modulo 
4, and take binary vertex states. Then y = (1/1,1/2,1/3,1/4), n[l] = (1,2,4), and y[l] = 
{Vii 2/2, Hi)- If we use the Boolean nor-function nor 3 : F2 — > F 2 (i.e. ECA # 1) given by 
nor 3 (a;, y, z) = (1 + x)(l + y)(l + z) to induce the vertex functions we get, e.g. Fi(y) = 
(nor 3 (?/[l]), ?/ 2 , 2/3, 2/4)- With update order 7r = (1,2,3,4) we get the SDS map 

(2.5) [Nor y , tt] = F 4 o F 3 o F 2 o F 1 . 

It is easy to verify that [Nory, vr](0, 0, 0, 0) = (1, 0, 1, 0). In contrast, for a parallel update 
scheme the state (0, 0, 0, 0) would have been mapped to (1, 1, 1, 1). The entire phase space 
of the SDS map in (12. 5p is displayed on the left in Figure [TJ 

What follows is a short overview of functional and dynamical equivalence of SDSs. The 
analysis is largely concerned with the aspect of update order and characterizes SDS maps 
for a fixed graph Y and fixed V-local functions 5V in terms of w. It will provide the basis 
for cycle equivalence. 



Also referred to as the Kleene star or Kleene closure of v[Y]. 
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2.1. Functional Equivalence. Two SDSs are functionally equivalent if their SDS maps 
are identical as functions. For a fixed sequence 5V, a natural question to ask is when is 
[3V,w] = [5V,?//] for w,w' G Wy. The update graph U{Y) provides an answer to this. 
The update graph of Y has vertex set Wy. Two words w ^ w' are adjacent if they have 
equal length, say m, and (i) they differ only by a transposition of entries k and k + 1, 
and (ii) {wk,Wk+i} e[Y}. The finite subgraph U(Y) of U(Y) induced by the vertex set 
Sy is called the permutation update graph, and is denoted U (Y) . Clearly, it is a union of 
connected components of U(Y). Both U(Y) and U{Y) are examples of neutral networks 
as mentioned in the introduction. The update graph [/(CirQ) is shown Figure [2j 

Let ~y be the equivalence relation on Sy defined by n ~y ir' iff 7r and n' belong to the 
same connected component in U(Y). We denote equivalence classes as [ir]y and the set of 
equivalence classes by Sy/~y, i.e. 

Sy/~y= {[ll]y | 7T G Sy} . 

By construction, ir ~y o implies the equality [^y,^] = [3V,c]. If the vertex functions 
are the Boolean functions as in Example 12.21 then [Nory,7r] = [Nory,a] implies tc ~y a 
(see |ll|]). It follows that |5y/~y| is a sharp upper bound for the number of functionally 
non-equivalent permutation SDS maps obtainable by varying the update order. 

Functional equivalence can also be characterized through acyclic orientations. An orienta- 
tion of Y is a map Oy \ e[Y] — > v[Y] x v[Y] that sends an edge {i,j} to either or 
(J, i). Let Acyc(y) denote the set of acyclic orientations of Y, that is, orientations that 
contain no directed cycles. In [l3j it is shown that there is a bijection 

(2.6) fy. S Y /~y — Acyc(r) . 

A permutation tt G Sy defines a linear order on v[Y] by tt^ = % < v j = iiz iff k < I. This 
order defines an acyclic orientation Y where O y ({v, v'}) equals (i>, v') if v < n v' and (v f , v) 
otherwise. The map fy in (I2.6P sends [k]y G S*y/~y to Y . It follows that 

(2.7) a(Y) = |Acyc(y)| 

is a sharp upper bound for the number of functionally non-equivalent permutation SDSs 
that can be obtained by varying the update order. The result can be extended to general 
word update orders w G Wy. We do not review this here, but refer to 



2.2. Dynamical Equivalence. Two finite dynamical systems 0, ip: K n — > K n are dy- 
namically equivalent if there exists a bijection h: K n — >■ K n such that 

(2.8) <j)oh = hoijj. 

With the discrete topology the concepts of dynamical equivalence and topological conjuga- 
tion coincide. Thus, the difference between functional and dynamical equivalence is that in 
the former, the phase spaces are identical, but in the latter, the phase spaces need just be 
isomorphic. Update orders that are related by an automorphism of the base graph give rise 
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to dynamically equivalent SDSs. The number of orbits ct{Y) under the action of Aut(F) 
on Sy/~y given by 7 • [tv]y = [7 * 7r]y, where 

(2.9) j*w = 7(101), . . . niwm) , 

is an upper bound for the number of SDS maps up to dynamical equivalence. This follows 
since for SDSs with Aut(y)-invariant vertex functions one has (see [Til ] ) 

(2.10) [3V,7*7i"] °7 = 7° [dv,^] 

for all 7r G Sy and all 7 e Aut(F). Via the bijection in (12.61) . this action carries over to an 
action on the set Acyc(Y), and the number of orbits is given by 

1 v ;i 7 eAut(r) 

Here (7) \ Y denotes the orbit graph of the cyclic group G = (7) and Y, see 0, 0]. This 
bound is known to be sharp for certain graph classes [3] , but in the general case this is still 
an open problem. 



3. Cycle Equivalence 



Definition 3.1. Two finite dynamical systems <\>: K n — ► K n and ip: K m — > K m are 
cycle equivalent if there exists a bijection h: Per(</>) — ► Per(-0) such that 

(3.1) V'IPer(^) h = h O 0| PerW , 

where i/' | p er (v) an d 0|p er (0) denote the restrictions of the maps to their respective sets of 
periodic points Per (ip) and Per(0). Two systems <fi and ip with identical periodic orbits are 
functionally cycle equivalent. 

Example 3.2. As an illustration we continue Example 12.21 with Y = G1T4 with vertex 
functions nor 3 : F| — ► F 2 with update orders n = (1,2,3,4), ir' = (1,4,2,3) and n" = 
(1,3,2,4). The two SDS maps [Nory,7r'] and [Nory,7r"] are cycle equivalent, which can 




Figure 1 . The phase spaces of Example 13.21 

be seen in the two rightmost phase spaces in Figure [TJ They are not functionally cycle 
equivalent. Later we show that for Y = G1X4 there are at most 2 cycle configurations when 
K = ¥ 2 = {0,1}. 



'In general one can define this where (j> and "0 have different vertex state sets K and K' . 
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It is clear that both functional equivalence and dynamical equivalence imply cycle equiva- 
lence. Define a, p G S m to be the permutations 

<x = (m,m-l,...,2,l), p= (l, m )(2,m-l)...([fl,LfJ + 1) , 

and let C m and D m be the groups 

(3.2) C m = (a) and D m = (a, p) . 

Both C m and D m act on the set of length-m update orders via (12.31) . Define the s-shift 
cr s (w) = a s ■ w, and the reflection p(w) = p ■ w = (w m ,w m -i, . . . ,W2,W\). We can now 
state one of the main results. 

Theorem 3.3. For any w G Wy, the SDS maps [$y, w] and [§y, c s (w)] are cycle equiva- 
lent. 



Proof. Set Pk = Per[Jy, crfc(w)]. By the definition of an SDS map, the following diagram 
commutes 

(3.3) P fc _ x P fc _, 



F-w(k) 



F-w(k) 



[SV, <r h (w)] 



Pk 



for all 1 < k < m — \w\. Thus we obtain the inclusion F w (k){Pk-i) C Pk, and since the 
restriction map F w ^y. Pk-i — ► F w ^)(Pk-i) is an injection, it follows that |Pjt_i| < \Pk\- 
We therefore obtain the sequence of inequalities 

\Per[$ Y ,w}\ < iPerfS^o-iWH < • < |Per[^ y , a- m ^(w)]\ < \Per[$ Y ,w]\ , 

from which it follows that all inequalities are, in fact, equalities. Since the graph and state 
space are finite all the restriction maps F w ^ in (I3.3P are bijections. Clearly (13. ip holds 
with h = F w (fy, and the proof follows. □ 

Theorem 13.31 shows that acting on the update order by the cyclic group C m preserves the 
cycle structure of the phase space. We point out that this result holds for any finite set 
K. For K = F 2 the cycle structure is also preserved under the action of D m , and is a 
consequence of: 

Proposition 3.4 ([U]). Let (Y, ^y,w) be an SDS over F 2 with periodic points P C Fg. 
Then 

(3-4) ([dv^w]]^- 1 = [d Y ,p(w)]\ p . 

This result follows since for each vertex function the restriction fi(—',y'[i])' F 2 — > F 2 
is a bijection for each fixed choice of y'[i}. There are only two such maps: the identity 
map yi i— > yi and the map yi \— > 1 + y^. From this it follows that composing the two maps 
in (13 .4p in either order gives the identity map, see [ill]. The next proposition is now clear: 
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Proposition 3.5. For K = F 2 the SDS maps [$y,w] and [$y,p(w)} are cyc/e equivalent. 

Thus, for any g <E G = C m the SDS maps [#y, w) and [3y,<7 • to] are cycle equivalent, 
where \w\ = m. If K — F 2 the same statement holds for C7 = D m . We now have the 
following situation: elements ir and n' with [tt]y 7^ [7r']y generally give rise to functionally 
non-equivalent SDS maps. If there exists g G G, n G [7r]y and ft' G [VJy such that g-jf = ir', 
then the classes \iz]y and [V]y give rise to cycle equivalent SDS maps. 

Let Star n be the graph with vertex set v[Star„] = {0,1,..., n} and edge set e[Star n ] = 
{{0,i} I 1 < % < n\. As a particular example we have: 

Corollary 3.6. Let Y = Star„ with n > 2. For a fixed sequence $y of Aut(Y) -invariant 
Y -local maps all permutation SDS maps of the form [$y, 7f] are cycle equivalent. Any SDS 
map of the form [Nory,7r] with ir G Sy has precisely one periodic orbit of size three, and 
2 n ~ l — 1 periodic orbits of size two. 

Proof. We have Aut(Star n ) = S n since the automorphisms of Star„ are precisely the ele- 
ments of Sy that fix the vertex 0. An orbit of Aut(Star n ) in Sy/^y contains all equivalence 
classes [ir]y for which the position of in tt coincide, i say. Thus for < i < n all permuta- 
tions that have vertex in the i th position give rise to dynamically equivalent SDS maps. 
Pick n = (0, 1, 2, ... , n). By Corollary I3.5j all permutations that are shifts of n give cycle 
equivalent SDS maps. The second par t now follows by inspection of one of the possible 
phase spaces. They are all listed in [llj], but without enumerations of periodic orbits. □ 

4. Combinatorial Constructions for Cycle Equivalence 

4.1. Neutral Networks. In the remainder of this paper we will only consider permutation 
update orders, although it is not hard to see how this can be extended to systems with 
general word update orders. To start, we define two graphs over S"y/~y whose connected 
components give rise to cycle equivalent SDSs for a fixed graph Y and a fixed sequence $y- 
Since cycle equivalence is a coarsening of functional equivalence, it is natural to construct 
these graphs using Sy/~ Y as vertex set rather than Sy- 

Let C{Y) and D(Y) be the graphs defined by 

V[C(Y)} = Sy/~y, e[C(Y)} = {{[n]y, [cT^jy} | TT G Sy} , 

v[D(Y)} = S Y /~y, e[D(Y)] = {{[n] Y , [p(n)] Y } | n e Sy} U e[C(Y)} . 

Define k(Y) and S(Y) to be the number of connected components of C(Y) and D(Y), 
respectively. It is clear that C(Y) is a subgraph of D(Y), and that 6(Y) < k(Y). By 
Theorem 13.31 is a general upper bound for the number of different SDS cycle equiv- 
alence classes obtainable through update order variations. For K — F 2 it follows from 
Proposition 13.41 that 8(Y) is an upper bound as well. 

Example 4.1. As in Example 12.21 let Y = G1X4. The permutation update graph f7 (Circ.4) 
has a(Circ4) = 14 connected components as shown in Figure [2j The graphs C(Circ4) and 
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1234 2341 1243 1423 1324 1342 

3412 4123 3241 3421 3124 3142 

1432 2143 2134 2314 2413 2431 

3214 4321 4132 4312 4213 4231 



Figure 2. The update graph [/(Circ 4 ). 

/}(Circ 4 ) are shown in Figure [3] where the dashed lines represent the edges in e[.D(Circ 4 )] \ 
e[C(Circ4)]. The vertices in Figure [3] are labeled by a permutation in the corresponding 




Figure 3. The graphs C(Circ 4 ) and -D(Circ 4 ). The dashed lines are edges 
in L>(Circ 4 ) but not in C(Circ 4 ). 

equivalence class in S"y/~y. The vertices of the cube-shaped component are all singletons in 
Sy/~y. The equivalence classes [1324]ci rC4 and [2413]circ 4 both consist of four permutations, 
while the remaining four vertices on that component are equivalence classes that contain 
precisely two permutations. Clearly, K(Circ 4 ) = 3 and <5(Circ 4 ) = 2. 

The following result gives insight into the how k- and 5-equivalent permutations are dis- 
tributed across the vertices of the update graph U(Y). 

Proposition 4.2. Let Y be a connected graph on n vertices and let g, g' G C n with g ^ g' . 
Then [g ■ 7r]y ^ \g' ■ n}y. If g, g' G D n with g ^ g' then [g ■ n]y = [g' ■ tt]y holds if and only 
ifY is bipartite. 

The proof, which can be found in [l(|, is by contradiction. We remark that if Y is bipartite 
if and only if \{[g ■ n]y : g G D n }\ = 2n — 1. 

4.2. Source-Sink Operations and Reflections of Acyclic Orientations. In this 
section we show how the component structure of C(Y) is precisely captured through 
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source-sink operations on acyclic orientations. The bijection in (12.61) identifies [%]y with 

Oy G Acyc(y). For any 7r G My, the orientation Oy 1 ^ is constructed from Y by 
converting vertex 7Ti from a source to a sink. Following [16( we call such a conversion a 
source-sink operation or a c/zcA;. It can be easily verified that this gives rise to an equiva- 
lence relation ~ K on Acyc(y). More precisely, two orientations Oy,0' y G Acyc(y) where 
Oy can be transformed into 0' Y by a sequence of clicks are said to be K-equivalent. This 
observation along with Theorem 13.31 shows that permutations that belong to ^-equivalent 
acyclic orientations induce cycle equivalent SDSs. By construction, the source-sink opera- 
tion precisely encodes adjacency in the graph C(Y), and the connected components are in 
1-1 correspondence with the ^-equivalence classes. Therefore, the number of equivalence 
classes in Acyc(F) under the source-sink relation equals k(Y), and is thus an upper bound 
for the number of cycle equivalent permutation SDS maps [3y, 7r] for a fixed sequence Jy- 

If K = ¥2 then Proposition 13.41 shows that reflections of update orders also induce cycle 
equivalent SDSs. On the level of acyclic orientations this corresponds to reversing all 
orientations. Through the bijection (I2.6P this identifies Y with the reverse orientation 
Y W \ the unique orientation that satisfies Y ({i,j}) 7^ Y n \{i,j}) for every {i,j} G 
e[y]. If two acyclic orientations are related by a sequence of source-sink operations and 
reflections, then they are said to be 5-equivalent. 

The notion of k- and ^-equivalence carries over naturally to update orders as well. Two 
update orders in Sy belonging to ~y classes on the same connected component in C(Y) 
[resp. D{Y)] are called ^-equivalent [resp. 5 -equivalent]. For two K-equivalent update orders 
7r and tt', there is a sequence of adjacent non-edge transpositions and cyclic shifts that map 
7r to 7r'. This is simply a consequence of the definition of Sy/~y and C(Y). 

We remark that from here there is a close connection to the structure of conjugacy classes 
of Coxeter elements, something we explain more in Section [7J The case of K = F 2 and 
reflections does not seem to play any role in Coxeter theory. 

5. Enumeration for k(Y) and S(Y) 

It is not difficult to show that S(Y) may be characterized in terms of k(Y). 

Proposition 5.1 ([lOj). Let Y be a connected graph. If Y is not bipartite then 5(Y) = 
\k(Y). IfY is bipartite then 5(Y) = \{k(Y) + 1). 

The proof uses the fact that p : Sy — > Sy extends to an involution 
(5.1) p* : Acyc(F)/~ K ^ Acyc(F)/~ K . 

The result now follows since p* has no fixed points if Y is not bipartite, and has precisely 
one fixed point if Y is bipartite. As a corollary, a connected graph is bipartite if and only 
if k(Y) is odd. In light of Proposition 15.11 we focus on the computation of k(Y) in the 
following. It can be shown that k(Y) does not depend on bridge edges, i.e., edges not 
contained in a cycle. 
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Proposition 5.2 ([9]). IfY is the disjoint union of graphs Y% and Y 2 , or if Y is a graph 
with e = {v,w} a bridge connecting the subgraphs Y\ and Y 2 , then 

(5.2) k(y) = «(yi)«(y 2 ) . 

For the computation of k(Y) we may therefore assume that Y is connected, and that every 
edge is a cycle-edge. Note that for the empty graph on n vertices E n we have n(E n ) = 1 
since a(E n ) = 1. The following corollary is immediate. 

Corollary 5.3. Let Y be a forest. Then k(Y) = 5(Y) = 1. 

^From Corollary 15.31 we get the following perhaps surprising results on dynamics of SDSs 
over trees: 

Proposition 5.4. Let Y be a forest and $v be a sequence of Y -local functions. Then all 
permutation SDS maps [3y,7r] are cycle equivalent. 

The proof is clear since ^-equivalent update orders induce cycle equivalent systems. So 
in other words, when Y is a forest, all permutation SDSs of the form [3y, tt] for fixed 3V 
share the same cycle configuration. This result may not be that significant if the SDS 
only has fixed points, or only has a small number of periodic points. However, for other 
functions, such as invertible ones, it is very powerful. The parity functions par fc : ¥ 2 — > ¥ 2 
are defined as par(y) = J2iV^ modulo 2, and are invertible for every graph Y (see fill]). 
Let Pary be the sequence of F-local functions induced by the parity vertex functions. 

Corollary 5.5. If Y is a forest then for any n,o- e Sy the maps [Pary,7r] and [Pary, a] 

are dynamically equivalent. 

The same result holds for the logical negation of the parity function, which is also invertible. 

We now give examples of the computation of k. Even though some of these results may 
be derived as special cases of more general results, they are needed for central examples 
in Section We begin with a result for k(Y © v), the vertex-join of Y and the vertex v. 
Recall that the graph Y © v is constructed from Y by adding to Y the vertex v and joining 
this new vertex to each vertex of Y. 

Proposition 5.6. If Y is a graph with e[Y] ^ 0, then 

(5.3) k(Y@v) = 25(Y ®v) = a(Y) . 

Proof. Each K-class of Acyc(Y © v) contains a unique acyclic orientation where v is a 
source [§]. It follows that there is a bijection between Acyc(F © f)/~ K and Acyc(F), 
hence (I5.3p . The complete graph on n vertices is simply the vertex-join of the complete 
graph on n — 1 vertices, and thus we get the following corollary. □ 

Corollary 5.7. Let K n denote the complete graph on n vertices. For n > 2 we have 
K{K n ) = {n-l)\. 
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Proof. There are 2^ 2 ) orientations of K n , and by the bijection in (12.61) . precisely a(K n ) of 
these are acyclic, and this is equal to the number of components of the update graph U(K n ). 
Since U(K n ) consists of the n! singleton vertices in Sy, ct(K n ) = n\. By Proposition 15.61 
K{K n )=a{K n _ 1 ) = {n-l)\. □ 



The quantity k(Y) is in fact a Tutte-Grothendieck invariant: 
Theorem 5.8 ([9]). Let e be a cycle-edge ofY. Then 
(5.4) k(Y) = k(Y:) + k(Y:) , 

where Y' e is the graph obtained from Y by deleting e, and Y" is the graph obtained from Y 
by contracting e. 

The proof of Theorem 15.81 is quite involved, and along with Proposition I5.2[ it implies 
that k(Y) = T(Y, 1,0), where T(Y,x,y) is the Tutte polynomial 13]. In contrast, it is 



well-known that the number of acyclic orientations of a graph satisfies a{Y) = T(Y, 2, 0). 

A graph Y has an n-handle if it is of the form Y = Y' U Circ„ where Y' and Circ n share 
precisely one edge. 




Figure 4. A graph Y with an n-handle. 



Proposition 5.9. Let Y be a graph with an n-handle where Y — Y' U Circ n . Then 
(5.5) k(Y) = (n - I)k(K') ■ 

Proof. Let e' = {v,v'} be the edge shared by Y' and Circ n and let e be the edge in Circ n 
incident with v. By applying Theorem 15.81 and Proposition 15.21 for bridge edges we obtain 

k(Y' U Circ n ) = k(Y') + k(Y' U Circ^) . 

Equation (15.51) follows through repeated applications of this process. □ 

As a simple, special case of Proposition 15.91 we obtain /t(Circ ra ) —n—1. Just take Y' to be 
the graph with vertex set vfY 7 ] = and edge set {{1, n}} in Proposition 15. 9[ 
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6. k(Y) as a Complexity Measure 

The number of possible orbit structures that one can obtain by varying the update order 
is a natural measure for system complexity. As we have shown, k(Y) is a general upper 
bound for this number, and so is 5(Y) in the case of binary states. Since these bounds 
are graph measures we can characterize complexity in terms of the GDS base graphs. As 
we have seen, bridge edges do not contribute to periodic orbit variability at all, and so it 
suffices to consider the cycles of the graph. As can be seen in the case of Circ n , increasing 
the size of a cycle does not contribute much, e.g. ft(Circ n+ i) = /t(Circ n ) + 1. However, from 
the result on graphs with handles it follows that even the addition of a minimal handle 
doubles the measure k, i.e. k(Y U G1X3) = 2n(Y), where Y and G1T3 share precisely one 
edge. The following example shows the effect on complexity that results from increasing 
the radius of the rules for elementary cellular automata. 

Example 6.1 (CA rule radius vs. periodic orbit complexity). We have seen that K(Circ n ) = 
n — 1. Thus, for any fixed sequence of radius-1 vertex functions the number of distinct 
periodic orbit configurations is 0(n). This changes dramatically for radius-2 rules. In this 
case the GDS base graph is Circ n2 with 

v[Circ„ >2 ] = {1,2, ...,n}, and e[Circ nj2 ] = {{i, j} I 1 < \i - j\ < 2} , 

with index arithmetic modulo n. The auxiliary graph C\rc! n2 is obtained from Circ n2 by 
deleting the edge {2, n}. The case n = 7 is illustrated in Figure El 



1 


1 


1 


7 2 












6 \ / 3 






5 4 


5 4 


5 4 



Figure 5. From left to right: The graphs Circ7, Circ7 j2 and Circ' 72 . 

For simplicity we set g n = K(Circ nj2 ) and c n = ft(Circ4 2 ). Successive uses of the recur- 
rence (15 .4p with edges e\ = {l,n} and e 2 = {l,n — 1} for both Circ„ i2 and C\rc' n2 gives 

c n = c n _i + 2c„_ 2 + 2™~ 2 , and g n = # n _ 2 + c n + 2c„_ 2 , 

where C5 = 18, cq = 46, g$ = 24, and g& = 64. These recurrence relations are straightfor- 
ward to solve with 

/t(Cir< 2 ) = [(3n - 5)2 n - 4(-l) n ] /18 and 

/?(Circ nj2 ) = [(2n - 6)2 n + 9 - (2n - 3)(-l) n ]/6 . 

Thus, by increasing the rule radius from 1 to 2 we see that the number of distinct periodic 
orbit configurations is 0(n- 2 n ). The corresponding bounds for 5 are easily obtained from 
Proposition 15.11 
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We have seen how non-trivial symmetries in the base graph give rise to dynamically equiv- 
alent SDS maps when the vertex functions are Aut(y)-invariant. Since dynamical equiv- 
alence implies cycle equivalence we can construct a bound k(Y) in the same manner as 
for a(Y). This bound k(Y) thus reflects the added cycle equivalence that are due to 
symmetries and that arise for Aut(F)-invariant vertex functions. 

We close with an example that illustrates this and the results of Theorem I5.8[ and Propo- 
sitions 15.21 and 15.61 and 15.91 

Example 6.2. Let Y = Q\ be the binary 3-cube, which has automorphism group isomorphic 
to S4 x Z 2 . It is shown in [3j that = 1862 and that = 54. Thus, there are 

at most 1862 functionally nonequivalent permutation SDSs over Q\ for a fixed sequence 
of vertex functions. Likewise, there are at most 54 dynamically nonequivalent Aut(Q|)- 
invariant permutation SDSs. It is known that the bound ot(Q\) is sharp, since it is realized 
for SDSs induced by, e.g. the nor 4 -function. 

The number of cycle equivalence classes is bounded above by and from the recursion 

relation (15.41) we get (with some foresight at each step) 



where Propositions 15.61 and 15.91 were used in the last step. Since Q\ is bipartite we also 
derive 8{Q\) = (133 + 1)/2 = 67, and thus in the case of K = F 2 there are at most 67 cycle 
classes for a fixed sequence of vertex functions. Straightforward (but somewhat lengthy) 
calculations show that k(QI) = 8{Q\) = 8. In conclusion, we have 



a{Q\) = 1862, a(Q 3 2 ) = 5A, k{Q\) = 133 , 8{Q\) = 67 , R(Ql) = 8(Q 3 2 ) = 8 . 



Thus if 3v is a sequence of Aut^D-invariant IMocal functions, then there are at most eight 
different periodic orbit configurations for permutation SDS maps [5y, tt] up to isomorphism. 
Moreover, because = 8{Q\) taking vertex states from K = ¥2 does not improve this 

bound. 

This example is only meant as an illustration, and a systematic treatment incorporating 
the analysis of the functions R and 8 for general graphs will be pursued elsewhere. 




$1) + <^ 



7. Summary 



In this paper we have shown how shifts and reflections of update orders give rise to se- 
quential dynamical systems with isomorphic periodic orbit configurations. We have also 
shown how to bound the number of periodic orbit configurations, and have derived several 
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properties of this bound k(Y). For binary states we have shown how S(Y) applies to give 
a sharper bound. Both quantities k and S serve as measures for dynamical complexity. 

We also have k(Y) = T(Y, 1,0) where T denotes the Tutte polynomial. There are other 
mathematical quantities counted by T(Y, 1,0), and thus by k(Y). For example, the set 
/Kcyc v (Y) consisting of all the acyclic orientation of Y with v as a unique source is also 
counted by k(Y), see p]. In fact, for each ^-equivalence class there is unique acyclic 
orientation with v as the only source. This allows one to construct a complete set of 
representatives for permutations realizing the possible periodic orbit configurations, see [||. 

The notion of source-to-sink conversions also shows up in the context of Coxeter theory 
(see, e.g. [g] for definitions). For a Coxeter group with Coxeter graph Y the number of 
conjugacy classes of Coxeter elements (see [iq ) is also bounded above by n(Y), e.g. fisj ] . 
In general it is not known if n(Y) is a sharp bound, but is is known for special classes of 
graphs such as Circ n as shown by Shi in [l6(. This connection between Coxeter theory and 
SDSs could potentially be very helpful in further exploring the properties of asynchronous 
GDSs. 

In this paper we have not explored the question of when k (and 5 when K = ¥2) is a 
sharp bound. That is, for an arbitrary graph Y does there exist a sequence of vertex 
function whose number of non-equivalent orbit configurations equals k(Y)7 Proving this 
would require one to construct such functions for any given graph. We have also omitted 
computational aspects related to cycle equivalence. Given two SDSs, what is the compu- 
tational complexity of determining if they are cycle equivalent? Related questions have 
been asked for, e.g. fixed point reachability in [2], but see also [1] for similar questions. 
Additional future work includes extending our results from permutations update orders to 
general word update orders as well as further exploring the effects of symmetries in the 
graph and the computation of the bounds R and 8 as illustrated in Example 16.21 
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